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Abstract 


In this paper the fracture problem of a composite plate 
which consists of perfectly bonded parallel load carrying lam- 
inates and buffer strips is considered. It is assumed that the 
fatigue cracks may appear and spread in main laminates or in 
buffer strips or in both perpendicular to the interfaces. The 
external load is applied to the plate parallel to the strips 
and away from the crack region. The problem is solved for fully 
imbedded cracks and for broken laminates or strips and the corre- 
sponding stress intensity factors are calculated. 


1 . INTRODUCTION 

In designing with high strength compos i te ; sheet materials 
the use of relatively low stiffness and; high toughness, buffer 
strips oriented parallel to the main 1 oad-carryi ng laminates has 
been under investigation for sometime.' The practical aim of this 
particular design practice is to improve the fatigue crack propa- 
gation and arrest characteristics of the structure. The fracture 
process may start as the initiation of a fatigue crack at a local 
imperfection in the load carrying laminate. The cyclic history 
of the applied loads permitting, it may. then be possible for the 
crack to enter the adjacent buffer strips as a result of subcrit- 
icaTor critical growth i.vr-F.ati gue cracks may also appear and prop- 
agate in the buffer strips. Knowing the' f atiguje and fracture 

.iV. 

characteristics of both materials, for studies relating to life 
and structural integrity it would therefore be necessary to have 

*This.work was supported by the National Aeronautics and Space 
Administration under the Grant NGR 3 9 - 00 7-01 1 and by the Nation- 
al Science Foundation under the Grant GK-42771 X . 
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a reliable analysis of the problem. Particularly in fatigue 
crack propagation studies it would be very helpful to have a 
technique for the calculation of corresponding stress inten- 
sity factors. 

Even though in actual structures the primary laminates 
and the buffer strips would be anisotropic, in this study, 
largely for reasons of analytical expediency, it will be as- 
sumed that both materials are isotropic and linearly elastic. 

The laminates and the buffer strips containing cracks of arbi- 
trary length perpendicular to the interfaces are assumed to be 
perfectly bonded and periodically arranged (Figure 1). The com- 
posite plate is loaded parallel to the interfaces away from the 
crack region. The case in which the propagating crack termi- 
nates at a. bimaterial interface is separately studied. For 
various crack geometries and for a certian material combination 
the stress intensity factors are calculated. In recent years 
the solutions of various special cases of this problem have 
appeared in literature. Referring to Figure 1, for example, 
the case of a single strip, (i.e.,y 2 = 0) 1S investigated in [1] 
and [2], the problem of two half planes bonded through a cracked 
strip -(i.e., h 2 = °°, b. = 0, a<h-j) is studied in [3] and [4], the 
solution to the same problem with the broken laminate (i.e., 
h 2 = “, b = 0, a = hj’) may be found in [5], The problem of two 
bonded half planes with the finite cracks terminating at or 
going through the interface was considered in [6] and [7]. 


2 



2. FORMULATION OF THE PROBLEM 


Consider the elastostatic problem for the two-dimensional 
composite medium shown in Figure 1. It is assumed that two sets 
of periodically arranged laminates or strips of finite Widths 
2h i and 2 h 2 are perfectly bonded along their straight boundaries. 
The strips 1 and 2 are assumed to contain symmetrically located 
cracks of lengths 2a and 2b, respectively and to have the re- 
spective elastic constants p-| , k-j and u 2 , k 2 where = 3 - 4v^. 
for plane strain and k. = (.3 - v.)/(l + v.) for generalized plane 

J J J 

stress, v., (j = 1, 2) being the Poisson's ratio. The external 

J 

loads are assumed to be such that the periodic nature of the 
problem is preserved and the crack plane is a plane of symmetry. 

The ‘sol ut ion- of th.e.actual. .tracti on-free crack problem may be 
obtained by means of the Usual superposition technique. In this 
paper we will therefore consider the singular part of the problem 
in which the self-equi 1 i brati ng crack surface tractions are the 
only external loads. Thus, referring to Figure 1, the solution 
of the problem may be obtained by solving the related field equa- 
tions for the two strips under the following boundary and con- 
ti nui ty condi ti ons : 

u j (0 ,y) = u 2 (0,y) , v T (0 ,y) . = v 2 (0‘,y), (0<y<~), ( 1 . a , b ) 

.■ : .._sv :st .a lxx ( 0 ? y-) = :. a 2 x.x ( a i. x .y(°>y) = a 2xy^ 0,y ^’ (°$y <0O )> (2.a,b) 
u l ( " h‘i *y ) = o, a ] xy (- h i»y) = o» (o < y < °°)» (3.a,b) 

u 2 (h 2 »y) = °» °2xy^ h 2 ,y ^ = ° 5 (° <y<0O )> ( 4 - a ,b ) 

°1 xy ^ x = ^ * ( - h-|<x<0), a 2xy ^ * (0<x<h 2 ), ( 5 . a , b ) 

< ^ 1 yy(x, 0 ) = - p -j ( x ) , ( -hi <x<a-h.| ) , v i ( x ,0) - 0, • 


° 2yy ( x ,0 ) = "P 2 ^ x ^ » (h 2 -b<x<h 2 ), v 2 (x,0) = 0, 

( 0<x<h 2 -b) , 

°kyy^ x ,00 ^ = °* a kxy ^ x = °» ( k=1 » 2 )» 

t ' • » ... 

■ * i , . \ .. • . 


(7 .a ,b) 
( 8 ) 


where, because of symmetry, only one quarter of a. strip from each 
set is considered. In (1-8) u^ and v^ are the x and y-components 
of the displacement vector in the strip k, (k = l,2), o k ^ ^ , (k = l,2, 
( i , j ) = (x,y)) are the usual stress components and p-j and p 2 are 
known functions. Defining the coordinates 


x 1 = x + h-| , x 2 = x - h 2 , y 1 = y = y 


2 * 


(9) 


the solution of the problem satisfying the necessary field equa- 
tions and the conditions (3), (4), (5), and (8) may be expressed as 

u.(x.,y) = - L-f Q co {l[f i ( s ) - — g i (s)Jsinh(x.s) 


+ x i g.j { s )cosh (x^s ) }eos (ys ) ds 


-1 / K i~ 1 


- f f 0 °° m.(r)r" l (- 2 ry )e“ ry si n(x. ,r)dr , 


i (x i ,y) = \ (s) + ~k — 9i (s)]cosh(x.s) 


K • +1 


+ x i g^ (s ) s i nh (x. s ) }s i n(y s )ds 
" - -l / lc i +1 


+ ^ / " m i (r)r"'(- 2 — + yr) e~ yr cos (x i r )dr , 


(i = 1 ,2) 


(10. a, b) 


where f . , g i , and m. * (i = 1,2) are unknown functions which may be 
determined from the remaining six boundary conditions given by (1), 
(2), (6), and (7). From (10) the stress components may be obtained 
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as 


1 

2u f 



2 °° 

°ixx^ x i ,y ^ = " f / o C ° [f i ( s ) cosh ( x i s > + sx i g i (s)slnh(x.s)]* . 

O , * ; \ 

*cos(ys)ds - - f Q °° m. (r) (1 -yr)e‘ yr cos(x l .r)dr, 

= |/ 0 “{[f i (s) + 2g • (s)]cosh(x-s) 


1 

2y i 


a ixy< x 1-*> = 


sx^g^(s)sinh(x.s)}cos(ys)ds 
I C m i ( r ) (l+yr)e" yr cos(x.r)dr , 
\ / 0 °°{ [f i ( s ) + g i (s)]s1nh(x i s) 

sx.j g.j (s)cosh(x..s)}sin(ys)ds 
f / 0 °° yrm i (r)e _yr sin(x.r)dn, 

(i = 1 ,2). 


(11 .a-c) 


If we now substitute from (10) and (11) into the four contin- 
uity conditions (1) and (2) and take the inverse transforms in the 
variable s, we obtain four algebraic equations for the functions 


f 1 > g| »- f£> and g 2 '. Sol vi ng these equations , f ^ and g^ , (i = 1,2) 
may be expressed in terms of m-j (r) and m 2 (r). The remaining two 
mixed boundary conditions would then give a system of dual integral 
equations for m-| and m 2 - Appendix A gives the express i ons for f. 
and g.j , (i = 1,2). Rather than deriving the system of dual integral 
equations, in this paper the mixed boundary conditions are directly 
reduced to a system of singular integral equations. Before doing 
this we observe that the conditions (6.b) and (7.b) are equivalent 
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to 


str 


°. (a.<Tx. |<h.), f_ a ] v i '(x i , 0 )dxi= 0 , 


(i = 1,2; a-! = a, a ? = b) 


( 12 . a ,b) 


If we now define 


G i ( x i ) = ^7 v v ( x i , 0 ) , (i = 1,2), 


from (12. a) and (lO.b) it is seen that 


a • 


m i ( r) = - ^—jqr f Q G i ( t ) s i n r t dt, (a-, 


(13) 


a , a q b , 


i = 1,2) 


(14) 


Thus, substituting from (14), (Al), (A2), and (11. b) into ( 6 . a) and 
(7. a) we obtain a pair of integral equations for., the new unknown 
functions G-j and G 2 . With (tk as defined in (14) for y++0 the sec- 
ond integral in ( 11 . b) becomes 

lim [- — / °° m-(r)(l + yr)e" ,yr cos (x . r)dr] 
y ->+0 % 0 


i / 1 


/n ' ( 


1 


TT ( K ^ + 1 ) 0 V t-X i t + X\. ' i 


)G • ( t ) dt , (i = 1,2) 


(15) 


Similarly, the functions defined by (A.l) may be expressed as 
F,,.(s) = ^ 7 +f z/ 1 ' [ M i j ( s , t ) - M. j (s,-t)]G.(t)dt, ■ 

(i = 1 , 2 ; j = 1 , . . ,4 ) (16) 


i J 
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where 


M i i(s,t> = - [ s ( h i - 1 ) + ~ 2 — -]e" s ^ h i -t ^, 

M i2 (s,t) = [ s ( h i - 1 ) - — ] e ” s ^ h i ~ ^ , 

M i3 (s,t) = s(h i -t)e" s ^ h i“ t ^ , 

M i4 (s,t) = [s(h.-t) - 1]e' s(l V t) , (i = 1,2): (17) 

If we substitute from (15), (16), (A2), and (11. b) into (6. a) 
and (7. a) and note that because of symmetry 

G i (x i ) = - G.(-x.j), (i = 1,2; - a i <x i <a i ) (18) 

we obtain 

_ G , ( t ) . h 

f -a t^x — dt + U ^ c 1 "l (x-jtjG-j (t‘)dt + / b D k 12 ( Xl ,t)G 2 (t)dt 

TT (, 1 +<-, ) 

= - ' - 4 - y ' p 1 (x 1 -h- 1 ), ( -a<x-| <a ) , 

b G 2 (t) a b 

-b t - ~ x ~ 2 ' ~ dt + -a k 21 ( x 2 > t ) G i( t ) dt + / b k 22 (x 2 ,t)G 2 (t)dt 

tt(1+k 2 ) . 

= - -j— P 2 ( x 2 + h 2 ), (-b<x 2 <b) (19, a, b) 

where the kernels k..(x.,t), ( i , j = 1,2) are given by 

. * J * 

k i j ( x i , t ) = / 0 “ K i j ( x i ,t,s)e" s ^ h i _t ^ds. (20) 

The expressions for the functions K.. may be found in Appendix B. 

^ J 

The index of the system of singular integral equations (19) is +1 ; 
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therefore the solution will contain two arbitrary constants [8] . 
Theoretically these constants are determined by using the single- 
valuedness conditions ( 1 2 . b ) . However, due to symmetry, if one 
only considers solutions satisfying (18), (12. b) will be auto- 
matically satisfied. For a <h^ and b < h 2 (Figure -1) from Appen- 
dix B it may be shown that the kernels k..(x.,t), (i,j = 1,2) 

* J ' 

are bounded and continuous in, the corresponding closed intervals. 
Thus, in this case, (19) is a system of ordinary singular in- 
tegral equations with fundamental functions 

W] (t) = (a 2 - t 2 )' H , w 2 ( t ) =• (b 2 - 1 2 )~ h (21. a, b) 

: I 

and may be solved in a straightforward manner. 

All the known special cases may be recovered from (19) by 
letting the length parameters a, h ^ , b, and h 2 go to proper limits. 
For example, for h^ -*°° and k^ . vanish, (19) uncouples giving 

the integral equations for center-notched infinite plates. For 
b = 0, l^-* 00 , (19) reduces to an integral equation in G-|(t), the 
Cauchy kernel and k^ being the only kernels. From Appendix B, 
the integrand giving k^ may be shown to reduce to 

K-j -j ( x.| ,t ,s ) = (l-4a-|h.|Se ^1 S ) ^e ^l s * 

*{a-| cosh ( x ,s ) [ 1 -(* 2 ( 3 + 2h ^ s ) e ' 2 ^ 1 s ] 

-2x-| sa-| c^si nh ( x ,s ) e -2 ^ 1 1 s + [1 -2s (h-j -t)* 

*[cosh(x,s)(3-2h 1 s-a 1 e' 2h l s ) + 2x^ sinhfx^s)]}. (22) 
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where the constants and a 2 are defined in Appendix B. 

This is the kernel found in [5]. If y 2 = 0 (18) further reduces 
to the kernel obtained for the center-notched strip in [2]. One 
may also observe that the kernels k, . have the expected property 

* vJ 

of i nterchangabi 1 i ty in the sense that k 2 j may be obtained from 
k.|j by changing the necessary indices and letting and 

a 2~* a 4 * 

3. THE LIMITING CASE: BROKEN LAMINATES . 

In the limiting case of broken • laminates,, ( i . e . , a = h -j or 
b = h 2 ) the integral equations (19) are still valid. However, in 

this case the kernels k.., (i,j = 1,2) will not all be bounded in 

J 1 1 

the corresponding closed intervals. For example if a = h^ , txh 2 , 
examining the integrands K..(x.,t,s) given' in Appendix B it may 

1 J 1 

be seen that the kernels k -j 2 , k^ , and k 22 are bounded for all 
values of x^ , x 2 , and t, but k^ becomes unbounded as x^ and t 
approach the end points +h^ s imu 1 taneously . The important conse- 
quence of this is that for a = h-j (T9) is no longer a system of 
ordinary singul ar i ntegral equations, it has generalized Cauchy 
kernels, the fundamental matrix of the system is different than 
that given by (21), and the solution requires more care. Noting 
that since K-j ^ is bounded and continuous everywhere in 0<s<°°, 
referring to (20) any divergence in the integral must be due to 
the behavior of the integrand for s->°°. Thus separating the 
asymptotic part K^-j $ of the integrand K-j -j we find 

Kn (Xi ,t,s) =.K ns (x 1 ,t,s) + K 1 1 ^ ( x 1 ,t ,s ) , 
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k n ( X] ,t) - k i -j $ ( x i ,t ) + k -j i ^- ( x i ,t) , 


k lls( x l ,t ) = / 0 °° l< lls ^ x 1 ,t,s ^ e " S ^ hl . t ^ ds 5 
k llf( x l ,t ) = / 0 °° K llf^ x l ,t,s )e" S h] ^^ds. 


-h n S 


K n s ( x l »t» s ) = e 1 {a-|COsh(x-|S) + a 2 [1 - 2s ( h-j - 1 ) ] [ ( 3 


2h^)scosh(x-|S) + 2sx-|Sinh(x^h)]} 


( - h ^ < ( x | , t ) < h -j ) , 


( 23 . a-e ) 


where is bounded in -h^ < (x-| ,t )<h^ . The singular part of 

the kernel may be evaluated by. using the fol 1 owi ng resul t [10] 


00 

'o 5 


r> -s(2h.-t) f s '' nh(s;( l> = _cT . » -s<2h,-t), sinh < 5x l> 


cosh ( s x ) 


d r- 1 


{, 


dt 
X 1 


n o 


cosh ( sx-j ) 


} ds 


dt n "(2h ] - t) 2 -x-, 2 2h l _t 


}] , (ni = 1,2,..), (24) 


giving 


k 1 1 s ( x i »t) ( oii 3a 9 ) 


2h^ -t 


2a, 


2 '(2h 1 -t) 2 -x ] 2 [ ( 2 h - 1 ) 2 - x n 2 ] 2 1 


[x, t 


- ( 2 h -j - 1 ) 2 ( 4 h r - 3 1 ) .+ .[ ( 2h 1 -.t ) 2 -x.-, 2 ] _1 t 


- t) [(2h 1 -t) 2 (2h 1 2 -h ] t-3x 1 2 ) + x ] 2 ( 6h ] 2 -3h 1 t-x ] 2 ) ] } 


2 

{2a 2 (h 1 + x-j ) 2 ^ - d ~ 2 + 6a 2 (h.j + x-j ) d 


d x 


dx 


1 


1 


a, - 3a 9 -i 

— 2 -~ } [ t - ( 2 h ] + x ] ) ] ” 1 


10 



+ {2a 2 ( h-| -x j ) 2 2 ~ 6ot2(h-|-x-|) ^ 


dx 


dx 


1 


1 


— 4 — - } [ t -(2h-j -x-j ) ] " T , *( r h 1 <(x 1 ,t)<h 1 ) 


( 25 ) 


Substituted into (19) and together with 1/(t-x), 
gives a typical generalized Cauchy kernel. In ; t h i s case the 
fundamental solution matrix of (19) is of the form 


W] (t) = (h^-t 2 )^, w 2 ( t ) = (b 2 -t 2 .)" B , (0<Re(y,3)<l ) . (26. a, b) 

’i 

The characteristic equations giving the exponents y and 3 may be 
obtained from (19) by letting 

G,(t) = H,(t)(h, 2 -t 2 )- Y , G 2 { t ) = H 2 (t)(b 2 -t 2 ) :B (27. a, b) 

and using the standard complex function technique [8, 2, 5], In 
(27) H-j and H 2 are assumed to be bounded in the intervals. 

-h-|<t<h-| and -b<t<b, respectively. .Omitting the detailed manipu- 
lations (see, for example, [2], [5-7]) the characteristic equa- 
tions are found to be 

2c0STry + 4ct2(y-l)^ ~ (a^ + o^) = 0, 


C0S7rB " °- (28. a, b) 

We note that (28,a)is the same as the characteristic equation 
found in [6] by two different methods for a crack tip terminating 

at a bimaterial interface and (28. b) gives theiknown result 

. .. ! 

6 = 0.5. Similar results may be obtained for b = h 2 » a<h^ . 

1 1 


4. THE SOLUTION AND THE RESULTS 


For a < h i , the system of integral equations (19) has 

the fundamental functions as given .by (21) which are the weights 
of Chebyshev polynomials. Hence the system may be solved in a 

\ r 

straightforward manner by using the . Gauss -Chebysliev integration 
technique described in [9]. For a = h-| , b<h£ (or a<h^ ,b = h^) , 
the fundamental functions of the integral equations are given 
by (26) where 6 = 0.5 and y is the root of (28. a) for which 
0<Re(y)<l. It may be shown that for all possible material com- 

V 

binations y is real. In this case' the fundamental functions 
are the weights of Jacobi polynomials and the solution. of the 
integral equations may be obtained by using the Gauss-Jacobi 
integration technique (see, for exampl e , (5-7 ] ) . 

After solving (19) from the formulation given in this 
paper it is clear that all the desired field quantities may be 
expressed in terms of definite integrals with appropri ate Green 1 s 
functions as kernels and G-j and G^ as density functions. In 
fracture studies one is usually interested only in the stress 
intensity factors which may be defined in terms of cleavage 

i 

» 

stresses °-jyy ( x -j *0) » (i = 1,2) and may be expressed in terms 
of the derivatives of the crack surface displacements G-j and G 2 

as foil ows [6] : 



1 2 



where 


4y, 

a<h-| : k a = 1im ^(x^a) a ]yy (x ls O) = - lim /2(a-x 1 } G^x-,), 

x^-a 2 x-^a 


a = h, : k = 1 i m /? (x 0 +h 0 ) Y a, uu (x, ,0)= -2y*lim /? (h,-x, ) Y . 

1 3 x^-h 2 2 2 . 2yy 2 x T h 1 1 1 

•G j ( X] ) , (29 . a-c) 


^2 ' 1 +2^ ( 1 -y) 1 - 2ot 2 ( 1 - Y ) 

simry ^tic-jUg + y 2 +< 2 vi 


(30) 



As an example we consider the composite medium shown in Fig- 
ure 1. It is assumed that the problem is one of "plane stress", 

there is no constraint in x-direction, and the , plate is loaded 

( *) 

in y-direction sufficiently far away from the crack region. ' 

'■ i 

Thus, in the perturbation problem the crack surface tractions 

I 

are constant and satisfy the following condition: 


(p 1 /p 2 ) = e:,/e 2 


(31) 


where E^ and E 2 are the Young's moduli. 

! 

Figures 2-5 show some of the .calculated results. The ma- 
terial combination used (v-j = 0.35, v 2 = 0.45, = 6.65y 2 ) is 

assumed to approximate boron epoxy sheets having buffer-strips 


(*)This means that the dimension of the plate in y direction is 
large compared to that in x- direction. If the opposite is 
true and if the external load is applied through fixed grips, 
then in the uncracked plate e = 0 and the 2 ^atio of th| crack 
surface tractions becomes (pi/p 2 )=[E 1 (1 -v 2 )]/[E 2 (1 -vi )]. 

In either case the analysis does not apply to the strips on 
the sides. 
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of the same material but different stiffness. Figure 2 shows 
the stress intensity factor'k for b = 0 and varying ratios 

a 

a/h^ and h 2 /h.|. Referring to Figure 1., for a.= h^ (and 
0 5 b < h 2 ) the asymptotic value of the cleavage stress around 
the singular point is (see (29. c)) ; 


°2yy ^ x 2 ^ 


(x 2 +h 2 ) 


■y » (x 2 - x - h 2 ) 


(31) 


where, for the materials under consideration y = 0.7015. Figure 
2 also shows the normalized k for a = h, , b = 0. As expected 

a I 

for h 2 *>0 k a (a = h^ )-><». Note that for h 2 - 0 the problem reduces 
to that of a homogeneous plate with periodic collinear cracks. 
Figure 3 shows the similar results for a = 0, and b/ h 2 and 
h ^ / h 2 variable. In this case for b = ,h 2 the stiffness of the 

adjacent medium is higher and consequently the power of the stress 
singularity y becomes less, than 0.5 (y = 0.3621). The results 
given in Figures 2 and 3 indicate that the stress intensity 

factor k g for a crack in the stiffer material is higher and k^ 

\ r : r . 1 

for a crack in the less stiff material is lower than the stress 
intensity factor in the corres pond i ng homogeneous plate. 

Figures 4 and 5 show the results for the case in which both 
sets of strips contain cracks. These examples pre given for the 
same material pair as used in Figures 2 and 3, and for a fixed 
width ratio h ^ / h 2 = 4. Figure 4 shows the stress intensity 
factor k a for various values of a/h^ and varying b / h 2 . Again 
for a = h-| the power of the stress singularity y is greater 
than 0.5 and the scale for k fl is given on the rp'ght with a 
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normalizing factor Figure 5 shows the stress intensity 

factor k b for the same problem. The scale for k b /p 2 /F for the 
limiting case of a = h^ is shown on the right. Note that for 

a = h l as ancl for a<h l k b tends t0 zero as b ^" h 2’ 

This latter result is due to the change in the power of stress 

singularity from 0.5 for b<h 2 to y= 0.3621 <0.5 for b = h 2> 
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, APPENDIX A 


The expressions for the functions f. and g ^ , (i = 1 , 2 ): 

' . i- . ' 

The functions f. and g.,. (i '= 1 , 2 ), may be expressed in terms 
of m-j and m2 by substituting from (10) and (11) into (1) and (2). 
F i r s t d e f i n e 


2s , 00 


2 . K i‘ 3 , 2. 2,, sin(h i r) 


F n (s) = ^ S~ n j (r)[ 2 s‘ + (r‘+s‘)J 


(r 2 + s 2 ) 2 


dr , 


o,2 „ o K.+1 2 ? cos(h.r) 

F i2 (s) = ~ f o m i (r) [2r + ~2 (r +S )] r > (r 2 +s 2 ) 2 dr5 


cos (h . r) 


F i 3 (s) " ~T~ f o rm i (r) (r 2 + s 2 } 2 dr ’ , ; 

, o s i n ( h . r ) 

F i 4 ( s ) = / 0 °° r 2 m.(r) dr ’ (1 = I 1 ’ 2 > * (Al.a-d) 


Equations ( 1 ) and ( 2 ) would then give 
4 

fl (S) j=l (a i j F l j + b i j F 2j ) » 
g i (s) - ? =1 ( c ij F lj + d ij F 2j ) * l* 


= 1 , 2 ) , 


( A2 . a ,b) 


where 


a ^ 1 ( s ) = [n-j cosh ( h 2 s ) - c^T^/S-j, 
a 12 (s) = [ n T sinh ( h 2 s ^ " c 12 T 1^ S 1* 
a 1 3 ( s ) = [-X^ngSi nh(h 2 s) - c i3 T -|l/ s i» 


a -j 4 ( s ) = [ A-| n 2 c °s h( h^s ) - c 14 T 1 ]/S ] 
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t> 1 1 ( S ) = a n (s) b 12 (s) = f-n 1 sinh(h 2 s) - d^T^/S^ 
b i 3 ( 5 ) = f n 2 s i nh ( h^s ) - d^T^/S^ 
b i ^ ( s ) = [n 2 cosh(h 2 s) - d^t^l/S^ 
a 21 (s) = .,[n 3 cosh(h 1 s) - c 2 1 T 2 i / s 2 , 

I 

a 2 2 ( s ) = [-n 3 sinh(h 1 s) - c 22 T 2 l/S 2 , 

a 23 ( s ) = [ n 4 s inh(h-jS) - c^Tgl/Sg, 

a 24 (s) = [n 4 cosh(h 1 s) - c 24 T 2 ]/ s 2 , 

b 21 (s) = a 21 (s), b 2 2 ( s ) = [n 3 sinh(h 1 s) - d' 22 T 2 ]/S 2 , 

b 2 3 ( s,) - [-X 2 n 4 s i nh( s ) - d 23 T 2 ]/S 2> 

b 24 (s),= [X 2 n 4 cosh(h 1 .s) - d 24 T 2 ]/S 2> (A. 3) 

where 

y 1 / y 2 ’ A 2 “ ^2^1 » 

X-jKp + 1 

5 1 = [(X 1 - 1 ) h 2 s - — 2 sinh(2h 2 s)]sinh(h 1 +h 2 )s , 

t l 

^ K 1 I 

5 2 = [(A 2 - 1 ) h i s - — 2 si nh(2h-j s ) ] si nh(h^ +h 2 ) s , 

i X i K-+1 

T i = t ( X ^ - 1 ) h 2 s - — 2 sinh (2h 2 s][h^s cosh(h^+h 2 )s 

+ %sinh(h 1 +h 2 )s] + [- — — — sinh(2h 2 s) 

- (X 1 +< 1 )h 2 s] Jgsi nh(h 2 -b 1 )s , : 

X « K -I + 1 

T 2 = [(X 2 “l)h^s - 2 — si nh(2h 1 s ) ] [h 2 s cosh ( h.| +h 2 ) s 


1 7 


I 



+ %si nh (h i + h 2 ) s] + hs i nh ( h-j -h 2 ) s [- ( A 2 +k 2 ) h-j s 

^ 1 ^ ~ fC 

1 + ~^~2~ “ si n h ( 2 h -j s) ] . (A. 4 ) 

c -j i ( s ) = L-|Cosh(h 2 s) + L 3 cosh(h-jS), 

i ■ 

C 1 2 ( s ) = L-|Sinh(h 2 s) - L ^ s inh(h-jS) , 

C 1 3 ( s ) = -L 1 sinh(h 2 s) + L3S i nh ( s ) , 
c 1 4 ( s ) = *--| cosh (h 2 s ) + A 1 L 3 cosh(h-| s ) , 

I 

d^Cs) = c -j -j ( s ) , d ] 2 (s) = - c -j 2 ' ( s ) j d -j 3 ( s ) = A 2 L ] s i nh ( h 2 s ) 

!' 

- L 3 s i nh (h-| s ) , 

i 

d -j 4 ( s ) = A 2 L -j cosh ( h 2 s ) + L 3 cos h ( h ] s ) , 
c 21 (s) = L 2 cosh ( h -j s ) + L 4 cosh(h 2 s), 

! 

c 2 2 ( s ) = -L 2 sinh(h-|S) + L^si nh ( h 2 s )., 
c 23 (s) = A ^ L 2 s i n h ( h 1 s ) - L^sinh^s), , 

c 24 (s) = A^L 2 cosh(h 1 s) + L 4 cosh(h 2 s), 1 

d 2 -| ( s ) = c 2 i (s) , d 2 2 ( s ) = " c 22 ^ s ^ d 2 3 ^ s ^ = " L 2 si r> h ( h i s) 

+ A^L^si nh( hgS ) 

d 2 4 ( s ) = t- 2 cosh(h|S ) + ' Agl^cos h ( h 2 s ) , j (A. 5 ) 

where 

1 K ? + l 1 ICt +1 

L 1 = p ~2 — s i nh (h 2 -h 1 ) s , L 2 = F ~i — sinhfhj -h 2 )s , 


18 



L 3 = ]J [ 0 -^2 > h 2 s ~ J 2(^2 +K 2) sinh (2h 2 s ) ] » 

! ' 

L 4 = 1 [(1 ~*i )h-, s ~ H( X i +k i )sinh(2h-|S)] , 

( K 1 + 1 ) ( Ko + 1 ) o ' 

D = - j s l n h ^ ( h i - h 2 ) s - C ( 1 - A -j 

K I A n 

- — 2 — ~ sinh(2h 2 s)][(7-A 1 )h 1 s - — ~ 


i 


i 




i 


i 

i 

i 

19 | 

i i 

» 

j 


) h 2 s 

sinh(2h-|S)] . (A. 6) 



APPENDIX B 


The expressions for the functions , (i»j - 1.2) which 
appear i n ( 20) : ^ 

-h, s 


K 1 1 (x 1 ,t ,s) 


'±L coshdx-js) [{(3-2h 1 s-(2h 1 s + 3)e _2s(h l +h 2 ) )Q ; 


Q 3 Q 8 Q9 

a^Q^Qg}{(2a2s(h-|-t)-a2 + e ^I^C^+O-ot-jO^Qg} 
Q g {(1+a 1 )Q 6 + "|+e~ 2h 2 s -2(h 1 -t)se" 2h 2 s }] 


+ e_h l S s-i'nh(x - 1 s ) tQ 3 (e _2h l S -a, 

i 

+ 2 0^2 s ( h i _ t ) } .+ ( 1 - ■ ot i ex 2 ) Q 8 1 » 


1 +a 


K 1 2 ( X "I ,t,s) " P 3 Q 8 Q g 


e‘ h l S ! cosh(x 1 s)'[{Q 3 (3-2h 1 s-(3 + 2h 1 s)e" 2s(h l +h 2 ) ) 


- 2 h„s 


■ a l Q l Q 5 HQ 5 (2a 4 s(h 2" t) ~ 2 )+Q 3 } 


- Q Q {4h 9 se' 4h 2 s +(l-e~ 4h 2 s )(2sh 2 -2st-l )}] 

i ! i 

1 +a 9 _K c t \ 

+ 2x-|Se"|l b S-inh(x 1 s){Q 5 [2a 4 s(h2-t)-a 4 

+ e ~ 2 h 2 s ] + Q-j) , • | 


K 21 (x 2 ,t.s)_» e- h 2 S |COsh(x 2 s)[{Q 4 (3-2h 2 s-(3 t 2h 2 s)e- 2s(h l +h 2 ) ) 


2h, s 


+ a 3 Q2Q 5 }{Q 5 (4 2 -2a2s(h r t)-e;^ n l )+Q 4 ) 

_ 1^1 Q{ 4hse _4h l S +(1-e _4h i S )(2h 1 s-2ts-l)}] 
I +01/1 y <1 


1 +a y 


2sx 2 e' h 2 S sinh(x 2 s) [Q 5 '( a 2 - 2a 2 s ( t ) -e~ c " 1 3 ) +Q 4 ] 


-2h, s 


22 (x 2’ t,s) = Q 4 Q ' q Q^ cosh(x 2 s) [{Q 4 (3-2h 2 ?-(3+2h 2 s)e 2s(h l +h 2 } ) 


-2h„s 


-2h,s 


- Qgd+e 2 -2s(h2‘t)e 1 +(1 +013 ) Q ? } ] 

2sx ? h c -?h <; 

+ q “ e" n 2 s sinh(x 2 s) [Q 4 (2a 4 s(h 2 -t)-a 4 +e ^ n 2 s ) 


- C -° t 3 0t 4 )Q 5 3 


where 


(B.la-d) 


cx 1 (K-|y 2 -< 2 yi )/(y2 +K 2 y l ) *\ ^2 

j 

a 3 = (< 2 y 1 -< 1 y 2 )/(y 1 +K 1 y 2 ) a 4 = ( y 1 -y 2 ) / ( y 2 + < 2 y 1 ) , 


4h 2 S -2h„s 


4h„s 


! 1 

( 3 


a 


Qo = e 


3 

-4h„s 


2 -1 +e 2 , ;Q 2 


4h l S e " 2h T s -l +e" 4h l s , 


a 


T 


Q 5 = e‘ 2h l s -e" 2h 2 s , Q 6 = 1 +4sh 2 e" 2s h 2" e , 2 , 
Q 7 = l+4h 1 se“ 2h l s -e" 4h l s > Q 8 = 1 -e‘ 2s ( h l +h 2 } , 


( B . 2 ) 


2 a + 4a 4 sh 2 e _2h 2 s -l ,! Q 4 = e _4h l s +4a 2 sh ] e' 2h l S -1 , 


Qg = ( “-j " 1 + ^ 3 ^ 4 * : 


( B . 3 ) 
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